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The integrable SM(4)-invariant spin-ladder model with boundary defect is studied using the Bethe 

ansatz method. The exact phase diagram for the ground state is given and the boundary quantum 

^ l- critical behavior is discussed. It consists of a gapped phase in which the rungs of the ladder form 

D ' singlet states and a gapless Luttinger liquid phase. It is found that in the gapped phase the boundary 

^AJ I bound state corresponds to an unscreened local moment, while in the Luttinger liquid phase the 

\^ . local moment is screened at low temperatures in analogy to the Kondo effect. 

75.10.Jm, 75.30.Kz, 75.40.Cx 
^ , I. INTRODUCTION 



Spin ladder systems are an active field of research in condensed matter physics experimentally realized in some 
quasi-one-dimensional materials .u The simplest example is the two- leg isotropic spin- 1/2 ladder which has a gapped 
'^ ,: ground state. Generalizations of ladders to more legs and couplings bej'OJid nearest neighbor exchanges show a 
G ■ remarkably rich behavior and can interpolate among a variety of systemsBcl For instance a dimerized phase dsiven 
O , by biquadratic interactions was predicted in Ref. [2] and then demonstrated in a generalized spin-ladder modeB by 
constructing the exact ground state and studying elRtnentary excitations. Recently, a modulation structure induced 
by frustration was reported in a two-leg spin ladder.B Quantum phase transitions from a gapped phase to a gapless 
phase were theoretically predicted -and experimentally studied in the Heisenberg ladder system Cu2(C5Hi2N2)2CIa 
*^ ' in the presence of a magnetic field .0 

L^ . Solutions of integrable models provide a useful starting point for the understandijic; of more general correlated many- 

^-^ body systems. A fewJntegrable spin-ladder models have recently been proposedJQ^EiJ For instance the integrable su{A) 



invariant spin laddeuJ|represents a special case of the Nesesyan-Tsvelikjnodel.El This model was recently generalized 

f^ ■ to the multi-leg casalil and hole-doping of the ladder was also studiedJla In general some four-spin interaction terms 

^^ ' must be included in integrable models (required by the integrability) , which could be either related to spin-phonon 

^P mediated interactions or in the hole-doped pha|e they could be generated by the Coulomb repulsion between the 

■ " holes moving in the spin correlated backgrounds The importance of a biquadratic interaction for properties of Cu02 

plaquettes was pointed out in Ref. [13] and its effect on excitations in a spin ladder was studied in Ref. [3]. Sora£ 

experiments indeed revealed that such multi-spin interactions are realized in-solid films of ^He absorbed on graphitc,Efl 

in a two-dimensional Wigner solid of electrons formed in a Si inversion layer Jlj in bcc solid '^iJeJlj and in heavy fermion 

^ ' systems. 

Q \ Impurities always play a relevant role in low-dimensional systems. ThisJs especially the case in Luttinger liquids, 

O ■ where an impurity may drive the system to a strong coupling fixed point |1I1 which corresponds to an open boundary 
^ ] condition at the impurity site for low energy excitations. Therefore, the boundary impurity is of pi aijtic ular interest in 
• ^ I quasi one-dimensional systems. Typical examples are the spin chain mth. boundary magnetic fields ,Pn^i or equivalently 
^\ two-dimensional classical statistical systems pith boundary fields,tj'L3 and a quantum impurity coupled to a one- 
dimensional strongly correlated electron host.OE3 

In this paper we study the open sM(4)-invariant spin ladder with a boundary defect. Boundary effects can arise 
for example if (i) the transverse coupling at the boundary is different from that in the bulk or if (ii) the rung-rung 
coupling at the boundary is different from that in the bulk. We will consider the following two model Hamiltonians: 
Model I 



N-l -, N 



il = - ^ (1 + (Tj • aj+i){l + fj ■ Tj+i) + jJ^^j ■ Tj + -J'ai ■ fi, (1) 



/I / J V J J \ -^ / \ ' J •' ' A / J J •'A 



Model II 



^N-1 ^ N 

H ^^^i^+^j- '?j + l)(l + -Tj ■ Tj+l) + ■^J'^^j ■ Tj 

+ ic/(l + ai-a2)(l + ri-T2), (2) 

where dj and tj are Pauli matrices acting on the site j of the upper leg and lower leg, respectively, and J represents 
the transverse rung coupling constant in the bulk, while J' is the coupling at the boundary rung. U denotes the 
rung-rung coupling strength between the first and second rungs. |Without the boundary defects (i.e., J' — J, U — 1), 
the model is exactly solvable with periodic boundary conditions .13 In Sect. II we show that model I is integrable for 
arbitrary J', obtain the corresponding Bethe ansatz solution, and discuss the physical consequences of the boundary 
bound state. In Sect. Ill the exact solution of model II is presented for arbitrary U. Conclusions follow in Sect. IV. 

II. SOLUTION OF MODEL I 
A. Bethe ansatz formulation 

The quantum states of a single rung are \a^,T^ >. It is convenient to define 

|o>=^(IT,i>-U,T>), 
|i>=|T,T>, |3>=U,i>, (3) 

|2>=^(IT,i>+U,T>), 

which satisfy the orthogonality relation < a\P >= 0. For simplicity we have omitted the rung index j and obviously, 
the first state denotes a singlet rung and the latter three the triplet states. We introduce the Hubbard operators 

X°'f^ = \a></3l a, /3 = 0,1, 2, 3, (4) 

and rewrite Hamiltonian (1) as 

N-l 3 N 11 

^ = E E ^f^'^i - ^E^f - J'^i" + \j(N - 1) + ij', (5) 

and the total number of a rungs can be expressed as N^ = X]i=i^i""- ^^ ^^^^ ^^Y ^^ have reduced model I to 
an su(4)-invariant spin chain with an effective magnetic field or, equivalently, to an sM(4)-invariant i — J mode£3a 
with a finite chemical potential and a boundary potential. Both the effective chemical potential J and the boundary 
field J' lift tbe-ipji^) symmetry of the Hamiltonian. Xiie Hamiltonian (5) can be diagonalized via either the algebraic 
Bethe ansatzE3E3l2l or the coordinate Bethe ansatz.Ea As the pseudo vacuum we chose the state in which all rungs 
are in a singlet, 

|f2 >= |0i >(g)|02 >®---® IOat > . (6) 

The Bethe wave functions can be constructed as 

where ^ is an amplitude, the sum over j„i runs from 1 to A^ and the sum over a™ from 1 to 3. The elimination of 
the "unwanted" terms from ll\^ > yields the standard nested Bethe ansatz equations (BAE): 

27V 




Ml _ \ _ 

TT TT /^" "^^2 2 TT /^" ' '^'^ 2 i^\ 



r=±7/5 



Vi -ri^y+i ^-^^ AX j^_j ^r^a, + ^ 



where the parameter 77 is determined by J and J' via 

J-J'-l^^, (9) 

1-77 

and Ml = Ni + N2 + N3, M2 = N2 + N3 and M3 = N3. Here Xj, jia and vs represent the rapidities of the flavor 
waves. The energy spectrum of the Hamiltonian (1) is given by 

A/i , „ 

^ = -E(t^-T-^)-4^^-4(^'-^) + ^-1- (10) 



B. Ground state properties 

For the periodic boundary conditions aU the rapidities are real in the ground state. For J = A the system has a 
quantum critical pointjj When J > 4, the ground state is the pseudo vacuum \Vl >, i.e. all rungs are in the singlet 
state. For J < 4, on the other hand, there is a continuum of excitations given by a Luttinger liquid of in general three 
components. 

The boundary defect may change the phase diagram close to the end point of the ladder. In particular, imaginary 
solutions of the BAE arising from the boundary scattering matrix (first factor on the right-hand side of the first set 
of Eqs.(8)) correspond to wave functions that fall off exponentially from the boundary. In fact, A = irj/2 is always a 
solution of the BAE in the thermodynamic limit N —^ 00 for 77 > and 77 7^ 1 (77 = 1 is a singular point of Eq.(9) 
which corresponds to J' = ±00). This imaginary mode represents the boundary bound state corresponding to a 
triplet rung. A careful analysis of the energy carried by the imaginary mode yields that the boundary bound state 
is not always stable (occupied) in the ground state. We limit ourselves to the situation of antiferromagnetic coupling 
(J > 0). From Eq.(lO) we have that the energy of the imaginary mode is 

eb = J--^. (11) 

1 — 77^ 

The boundary bound state is stable if ei, < 0. Otherwise the imaginary mode represents an excited state. We have to 
distinguish the cases J > 4 from J < 4. (i) For J > 4, there is a critical line given by 



ric-^l-j, (12) 

which separates the spin singlet rung ground state (0 < 77 < 77c and 77 > 1) from the spin triplet ground state 
(j]c < 77 < 1) at the boundary of the ladder. For < 77 < 77c and 77 > 1 we have et > 0, while for rjc < r] < 1 eb 
is negative and the boundary bound state is filled, (ii) For < J < 4 the bulk corresponds to a three component 
Luttinger liquid with real A, fi and v modes. Here the boundary bound state is stable in the whole parameter region 
< 77 < 1, but it is empty if 77 > 1. The boundary triplet state for < 77 < 1 is coupled to the continuum giving rise 
to a Kondo-like screening. 

Consider now the response of the system to an external magnetic field h. The magnetic field couples to the ladder 
via the Zeeman effect, i.e. the Hamiltonian has an extra term —hj^j^ii^j^^^j^)- The critical line J = 4 separating 
the gapped and gapless regions is now shifted to J — h — A. The bound st ate energy is also reduced by h. Hence, 
for J — h> A, the critical line Eq.(12) is now given by rjc — \J\ — 4/(J — h). The threefold degeneracy of the triplet 
rung state at the boundary is lifted by the field, so that an arbitrarily small field induces a finite magnetization ±1 
(depending on the direction of the field) due to the stabilization of the bound state. Therefore, the boundary quantum 
phase transition at r/ = 77c is of first order and the susceptibility is divergent at T = 0, following a Curie law. 

To summarize, the boundary phase diagram for J > ft, is shown in Fig. 1(a) and consists of seven regions. For 
J — h > A the bulk is gapped and we have argued that for 77 < there is no bound state and the rungs all form 
singlets. For < 77 < 77c the bound state is not stable and all rungs are in singlet states. For 77c < 77 < 1 there is a 
triplet state at the boundary (the triplet wavefunction falls off exponentially into the bulk of the ladder), and finally 
for 77 > 1 the bound state is again unstable. For J — h < A the bulk is a Luttinger liquid, without a bound state for 
77 < 0, with a stable bound state for < 77 < 1 and with an unstable bound state for 77 > 1. Below we show that a 
Kondo-like screening occurs for < 77 < 1. 

It is also interesting to study the situation for h — J > 0. The bulk is then always a Luttinger liquid and only 
four cases for the boundary bound state have to be distinguished (see Fig. 1(b)). We assume here that ft, — J is 



sufRciently small so that the spin ladder is not spin polarized. If 77 < there is no bound state, for < 77 < 1 the 
bound state is filled with a predominantly spin-up triplet state (the Kondo screening is quenched by the magnetic 
field), for 1 < 77 < ?7c there is an empty bound state, and for rj > rjc the bound state is again stable (with magnetic 
field quenched Kondo screening). 

The thermodynamics of the boundary defect can also be derived from the BAE, Eq.(8), following the standard 
method.L^tZI The thermodynamic BAE allow us to study the boundary quautum critical behavior. The boundary 
defect induces a "ghost spin" 77. However, unlike in su(2)-invariant models,E3 the ghost spin does not lead to an 
anomalous remnant entropy because the su(4) symmetry in the present model is already lifted by the finite J . 

C. The quantum critical line J = A 

Along the quantum critical line J = 4, the boundary defect can show critical behavior as the bulk doesE We now 
consider the case 77 > 1. In zero magnetic field, the ground state of the bulk consists only of singlet rungs. In a weak 
magnetic field some triplet rungs with S^ — \ appear in the ground state, while N2 and N^ still remain equal to zero, 
since the energy of the state |2 > is unchanged and that of the state |3 > is increased {h > 0). We denote with p{X) 
the distribution of real A modes, including the boundary-defect contribution. ^From the BAE, Eq.(12), we obtain 

pW+I dA'a2(A-A')p(A')=ai(A)-^a,(A), (13) 

where a„(A) = n/2TT[X'^ + {n/2)'^] and A^ = 1/(4 - h) - 1/4. For h « 1, we have A « \/7i/4 and Eq.(13) can be 
solved by iteration. 



TT Trr/N J \ TT 



with the ground state energy given by 



E = J^^dx(^4- ^^ - h^ piX) -^JN -^{J' - J) + N - I . (15) 

Combining Eqs.(14) and (15) we obtain the susceptibility of the system 

d^E [2 1 \ /l, 1 1 A ^,, 1, 

X = -ITTT = - - -777 ih-^ -^]+ 0{h^)- (16) 



dh'^ \TT irrjN J \4 37r 

The susceptibility diverges with the square root of the field as a consequence of the van Hove singularity of the empty 
A band. The boundary bound state removes one degree of freedom from the bulk, so that its contribution to the 
susceptibility is negative. This result is not surprising, since in this case J' is much larger than J. Consequently the 
boundary rung is in a tight singlet and hence insensitive to the field, so that the whole susceptibility is reduced. 

Similar arguments indeed yield a positive boundary susceptibility for 77 < 0. As discussed before, in the region 
< 77 < 1, a stable boundary bound state occurs and a small field already induces a finite magnetization. 

With a simple scaling approach we find that the low temperature specific heat and the magnetic susceptibility of 
the boundary defect at the line J = 4 behaves as 

<5C(T)^r^ <5x(r)^r-i (17) 

Such a result can also be predicted by a simple spin wave theory with a dispncsion relation e(fc) ^ fc^ or alternatively 
exactly via the low temperature expansion of lie thermodynamic BAEH3'LJ The boundary critical exponents in 
Eqs. (16-17) are exactly the same as for the bulk.Q 

D. Kondo effect in the gap less phase 

In the gapless phase, < J < 4, the system is a three component Luttinger liquid. In the sense of the su{A) t — J 
model, the triplet rungs are considered spin-1 hard-core bosons. The A rapidities represent the charge sector, while 
the p and v rapidities parametrize the spin degrees of freedom, which have sm(3) invariance. 



As discussed above (see Fig. 1(a)) a stable boundary bound state only exists at low temperatures for < 77 < 1. 
The boundary bound state corresponds to a local moment with spin 1. The boundary coupling J' does not break 
the SM(3)-invariance of the hard-core bosons and the boundary local moment is spin compensated in analogy to the 
Kondo effect. 

To show this we explicitly consider the imaginary mode 177/2 in the BAE, which then become 

A, -^ y ^ A, - §77 A, - z{l - f ) A, ~ z(l + 2) 
l^A, + 1 j A, + §77 A, +z(l - 2) A, +z(l + ^) 

Mi-l . . . Af 2 X , i 

X TT TT ^ -^^i-^ TT Aj - r/ia + 2 
iii^ A,-rA.+7llA,-.M.-r 

T-r TT- fJ-a - r^lp - I ^ -|-|- /ig - ^(l+r?7) -pj- /Jg - r-Aj - ^ -pj- Ma - '^'^a - 2 Qg^ 

ii/3i^"-'^^/5 + * ^4^*0 + 5(1 + ^'?) /=i Ma-rA, + § AiAig-rz., + §' 

n n :h^ - n n ^^^^^^^ 

r— ± 77^(5 



7^5 - rzy^ + 7 ^J-^ ^J-^ US - r^g + 



For large N the solutions of the BAE, Eq. (18), are strings of arbitrary length for all three sets of rapidities. We 
introduce the usual densities of A, /i, 1/ strings pi,„(A), p2,n{fJ'), Pz,n{v): and their respective hole densities p5'ti(A), 
P2 n(M)j P3,n(^)- In the thermodynamic limit these densities satisfy the following integral equationsLJ 

Pl,n(^) + zJ ^rnnPl,m{>^) = ^J BrnnP2.,mW + an(A), 

771 771 

PinW + Yl ^rnnP2.miX) = Y^ ^77777 (pi,m (A) + P3,77i(A)), (19) 

P3,nW + Y ^m«P3, 777(A) = ^ ^77777^2, 777 (A), 

777 7n 

where we neglected the boundary driving terms which are of order N^^. Here 

^77777 = [771 + 7^] + 2[77T. + 77 - 2] H h 2[|771 - Tl] + 2] + [|771 - 7l|] , 

7nin{m.7i} 

Bmn= Y Vn + n-2l + \\, 
1=1 

and [n] is the integral operator with kernel a„(A) and ao(A) is the (5-function. 
The free energy functional is given by 

^/^ = E [d\W,nPr,nW " T[p.,„(A) + p.^^JA)] ln[p.,„(A) + p,\(A)] 

r,77 "^ 

+ rp,,„(A) lnp,,„(A) + Tp'^.^iX) Inp'^JX)} , (20) 

where ei „ = — 27ra„(A) + n{J — /i), 62, „ f= .fj.n = Ti/i. Minimizing Eq. (20) with respect to the densities and taking 
into account the relations (19), we obtainHO 

ln(l + r;,,„) = ^+5]A-„ln(l + 77-i), r,s = l,2,3 (21) 

777, S 

with A"„ = AmnSr,s ~ Bynn(3T,s^\ + ^r,s-i), 77^,77 = Pr.nl Pr,n, and 77(^,\ = Tj^^^ = 0. An equivalent set of integral 
equations is 

ln77^,„ == G • [ln(l + rir,n+i) + ln(l + 77,-,,7-i)] - G * [ln(l + r;,7+i_„) + ln(l + 77;l\_„)] , 

ln??,M = - Y^(^)^'^.i + G*ln(l + 77^,2) - G* [ln(l + 77-_^\ i) + ln(l + 77-_\_i)] , (22) 

ln77i_„ J-h ln772.n ^. ln7;3.„ h 
hm = , hm = hm = — = 2a;o , 

77—^00 77 i 77— >oo 71 77— J-oo 77, J 



where • denotes convolution and G{X) ~ [2cosh(7rA)] ^. The equiHbrium free energy is 



F/N = -r^[n]ln(l 



Vil) 



At low T the exch ange J give s rise to a Fermi surface for the charges, which are only significantly populated in the 
interval |A| < A = \/l/J — 1/4, but are unoccupied for |A| > A. The low energy spin excitations, on the other hand, 
take place at very large rapidities (for h = the spin Fermi surface is at cxd). Hence, at low T the charge and spin 
sectors are well separated and only weakly coupled. Assuming complete decoupling of the spin and charge sectors, a 
solution of Eq. (22) can be easily obtained for large n 

sinh(na;o)sinh(n+ 1)2:0 _, ,„„, 

V2,n = ??3,n = r-T . , ,^ ^ 1 ■ (23) 

sinna;o sinn(2.To) 

Although the boundary bound state has also some charge fluctuations (triplet-singlet rung admixture), these do 
not affect the dynamics of the spin (weak coupling of spin and charge sectors) . We can then limit ourselves to the 
spin degrees of freedom of the bound state. Its contribution to the free energy is 

Fbs = -^TJ2 rfAK(A - \) + a„(A + \)] ln[l + ry^-i] (24) 

n=l-' 

or, equivalently, 

2 

Fts = FS, -\tY.Y.I ^^^^/(^ + ^'^^Z^) H^ + '?9+i.i] ' (25) 

r=± q=l •' 

where F^^ the ground state energy of the local moment and 

sin[7r(l - g/3)] 



G,{\) 



cosh(27rA/3) + cos[7r(l - g/3)] 



Substituting Eq.(23) into Eq.(25) we readily obtain that the residual entropy of the magnetic moment is exactly zero, 
which means that the boundary bound state is fully screened by tho-Kondo effect. These results are analogous to 
those of a Coqblin-Schrieffer impurity embedded into the same host.H 



III. SOLUTION OF MODEL II 



To show that model II is integrable we rewrite the Hamiltonian (2) as 

H — Ho + Hi, 

N-l 

Ho=Y.Pjj+i+UPi2, (26) 



J=2 

N 



H^ = -jt^r+T 



4 

where Pij = (1 + (7^ • (7j)(l -t-T^ • 7j)/4 is the permutation operator between rung i and rung j. Obviously, [H^, Hi] = 0, 
which means that they can be diagonalized simultaneously. We define the Lax operators 

S,,{\) = \-in, (27) 

which satisfy the Yang-Baxter relationa^J 

5y (A - ^l)SikiX)SJk{^J') = SJk{^J■)Sik{x)s^j{x - n). (28) 

As shown in Refs. [18,22], the monodromy matrix 

6 



TriX) = SNrWSN-lrW ' ' ' 52r(A)5l.(A - ^c)5l,(A + ic)S2r{X) ' ' ' Sn-It^SnAX) (29) 

satisfies the reflection Yang-Baxter equation 

Srr'iX - fl)Tr{X)Srr'iX + ti)Tr'{li) = T,,(/i)5,v(A + ^,)TriX)Srr' (X - fl) , (30) 

where c is an arbitrary constant and r and r' are indices of the 4-diniensional auxiUary space. From Eq.(30) we have 

[t{X),t{p)]^0 (31) 

with t{X) = trrTr{X). Therefore, t(A) serves as the generating functional of a variety of conserved quantities, i^o is 
related to t{X) as 

provided that U — 1/(1 — c^). Hence, Hq and Hamiltonian (2) are integrable in the sense of the algebraic Bethe 
ansatz. Following the standard procedure we readily obtain the BAE of model II 

\ ^j ~ l\ ^0 + *("^ + 2 ) '^J ~ ^\'^~ 2' TT TT Aj — 7-A; — 2 -p-|- Aj — r^a + 



Knn^^^^^n 



Aj- + ly Aj -i(c+ i) Aj +t(c- i) ;Lif_,-^ ^J -^^' +* ail Aj -r/x„ - f ' 

M2 . Ml \ i A/3 i 

-TT TT M« -r^lp-l ^ TT TT A^g - ^-Aj - 2 TT Ma " ^-^^^ - 2 /ggN 



A/3 . A/j 

J/5 — ri^-y — i TT TT '^'^ ^ ^/''" 



nn ' =nn 

r=±'^:^5 ' r=± a=l " ^" 



2 



with the energy eigenvalues given by 

A/i 

^ = -E(Ty-T-'^)-7'^^ + f^ + ^-2- (34) 

The Hamiltonian (32) is only Hermitian if the parameter c is either real or imaginary. For imaginary c, C/ < 1, and 
the first rung is weakly coupled to the bulk. For real c (we can consider c > because ±c yields the same U) the 
imaginary mode A;, = i(c — 1/2) is a solution of the BAE, Eq.(33), in the thermodynamic limit N ^ oo for c > 1/2. 
This bound state again corresponds to a triplet boundary bound state, i.e. its wave function falls off exponentially 
with the distance from the boundary. The energy of the bound state is e^ = J — l/(c — c^), which can be stable only 
if efc < 0, i.e. in the region 1/2 < c < 1. This corresponds to U > 1 for which the first rung is strongly coupled to the 
bulk. For c > 1, i.e. U < 0, the bound state has positive energy and is empty. The first rung is then fcrromagnctically 
coupled to the ladder. 

As in the case of model I, there is a critical line for J > 4 given by 

Jc = 7 2\ ' (35) 

(c — c^) 

which separates the region of an occupied (stable) and empty (unstable) bound state. Hence, when J > Jc the ground 
state is a spin singlet, while for 4 < J < Jc and 1/2 < c < 1 the boundary bound state is stable in the ground state. 
In summary, the physical properties of the two models are very similar. The main difFerenc£-|Can be understood in 
terms of the number of "ghost spin" solutions (number of bound state solutions) of the BAEJH3 which correspond to 
images of the real local moment. In the Kondo regime for J < 4 there are two "ghost spins" , c + 1/2 and c — 1/2, in 
the spin sector of the second model, while there is only one in model I. This can be read off from the impurity factors 
in the BAEs (8) and (33). Since the effects of the ghost spin contributions are additive, the physics in both situations 
is very similar. 



IV. CONCLUSIONS 

We studied two models for boundary defects of the open two-leg su(4)-invariant spin ladder. In model I the 
transverse coupling at the boundary rung is different from the bulk, while in model II the coupling of the first rung 
to the ladder is different from the rung-rung coupling in the bulk. The two models under consideration are integrable 
and we obtained the exact solution by means of Bethe's ansatz. Depending on the model parameters three situations 
may arise in the thermodynamic limit: (i) there is no imaginary mode solution of the BAE and hence the states of 
the first rung are part of the continuum of the bulk, (ii) an imaginary boundary mode exists, but corresponds to a 
positive energy et, i.e. the state is empty, and (iii) an imaginary boundary mode with negative energy exists. In 
case (ii) the boundary bound state does not affect the ground state properties, but does contribute to the finite T 
thermodynamics. The situation (iii) is the most interesting one, since a spin-1 boundary bound state is filled in the 
ground state. The boundary phase diagram of model I is shown in Figs. 1(a) and 1(b). 

The su(4)-invariant two- leg ladder has a critical line at J = 4. For J > 4 all the rungs of the bulk are in the singlet 
state, while for J < 4 the system is a Luttinger liquid of in general three components. For J > 4 a stable boundary 
bound state carries a magnetic moment (triplet state with wave function that falls off exponentially into the bulk), 
while in all other cases the boundary rung is in its singlet state. For J < 4, on the other hand, a stable boundary 
bound state carries a Kondo compensated (screened by the spin degrees of freedom in the Luttinger liquid) magnetic 
moment of spin-1, i.e. ultimately the ground state is a singlet. An unstable bound state just removes one degree of 
freedom from the Luttinger liquid. Both models considered here display similar properties at the boundary. 
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DMR98-01751 and No. DE-FG02-98ER45797. Y. Wang is also supported by the National Science Foundation of 
China. 
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Figure caption 

Fig. 1 : Boundary phase diagram of the two-leg su(4)-invariant spin ladder with different transverse coupling at the 
first rung (model 1) for (a) J — h > and (b) h — J > 0. The interaction strength at the first rung is parametrized 
by rj defined in Eq.(9). Quantum critical behavior with mean-field exponents is obtained along the line J ~ h = A. 
Several phases are possible at the boundary. In S all rungs are in a singlet state and there is no boundary bound 
state. In SI the ground state consists of singlet rungs but there is an empty boundary bound state. TBS refers to a 
phase in which the boundary bound state is stable, i.e. there is a boundary triplet state with wave function falling 
off exponentially into the bulk. LLl, LL2 and LL3 refer to a Luttinger liquid with no boundary bound state, with 
a stable boundary bound state (triplet state) and with an unstable boundary bound state, respectively. In (b) we 
assumed that h — J is sufficiently small so that the ladder is not spin polarized. 
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